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Power Divider 

Yu Christine Chen, Member, IEEE and Sairaj V. Dhople, Member, IEEE 


Abstract —This paper derives anaiyticai ciosed-form expres¬ 
sions that uncover the contrihutions of nodai active- and reactive- 
power injections to the active- and reactive-power flows on 
transmission iines in an AC eiectricai network. Paying due 
homage to current- and voitage-divider iaws that are simiiar in 
spirit, we baptize these as the power divider iaws. Derived from 
a circuit-theoretic examination of AC power-flow expressions, 
the constitution of the power divider iaws reflects the topoiogy 
and voitage profiie of the network. We demonstrate the utiiity 
of the power divider iaws to the anaiysis of power networks 
by highiighting appiications to transmission-network aiiocation, 
transmission-toss aiiocation, and identifying feasibte injections 
white respecting tine active-power flow set points. 

Index Terms —Power flow, transmission-ioss aiiocation, 
transmission-network aiiocation, feasibte injections. 


I. Introduction 

T his paper presents analytical closed-form expressions 
that map nodal active- and reactive-power injections to 
active- and reactive-power flows on transmission lines in a 
power network operating in sinusoidal steady state (see Fig. [T] 
for an illustration). We term these the power divider laws, 
since their form and function are analogous to and reminiscent 
of voltage- and current-divider laws that are widely used in 
circuit analysis ifTlFI 

Arguably, the most obvious application of this work would 
be in transmission-network allocation, where one seeks to 
equitably and systematically apportion the cost of loading 
transmission lines in a power system to constituent generators 
and loads. Indeed, we demonstrate how the power divider laws 
can be leveraged for this task. More generally, quantification 
of the impact of injections on flows is germane to numerous 
applications in power-system analysis, operation, and control. 
(We interchangeably refer to nodal complex-power injections 
and line complex-power flows simply as injections and flows, 
respectively. Furthermore, we implicitly consider injections 
to be positive for generation and negative for load.) In this 
paper, we provide a snapshot of additional possibilities by 
focusing on the particular tasks of allocating transmission¬ 
line losses and identifying feasible injections under the con¬ 
straint of line active-power flow set points. Through this 
promising (but by no means exhaustive) set of examples, 
we demonstrate how the power divider laws yield analytical 
insights, provide computational benefits, and improve accuracy 
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Fig. 1; The power divider laws: mapping injections to flows. 

compared to state-of-the-art approaches. For instance, with 
regard to the task on identifying feasible injections, using 
the power divider laws we formulate a linearly constrained 
least-squares problem, the solution of which returns the set of 
injections that respect line active-power flow set points. On the 
other hand, with regard to loss allocation, we recover exact 
expressions that map transmission-line losses to bus active- 
and reactive-power injections. This has been recognized to 
be a challenging problem in the literature 0, 0. Consider, 
e.g., the following remark from Q:"/---/ system transmission 
losses are a nonseparable, nonlinear function of the bus power 
injections which makes it impossible to divide the system losses 
into the sum of terms, each one uniquely attributable to a 
generation or load.” In this paper, using the power divider 
laws, we demonstrate how losses on an individual transmission 
line can be uniquely attributed to each generator and load. 

Given the obvious utility of mapping injections to flows in 
a variety of applications pertaining to power system operation 
and control, many algorithmic approaches and approximations 
have been applied to tackle this problem in the literature. 
Some relevant prior art in this regard include: numerical 
approaches a, Q, integration-based methods 11, utilizing 
current flows as proxies for power flows ||9l, and leveraging 
generation shift distribution factors Qol, d. Distinct from 
the numerical methods in 0-l9l, the power divider laws offer 
an analytical approach to precisely quantify the relationship 
between line flows and nodal injections. Furthermore, distinct 
from d, ini, our approach tackles the nonlinear power- 
flow expressions and acknowledges injections from all buses 
in the network; no assumption is made about the existence or 
location of a slack bus that makes up for power imbalances. 

The derivation of the power divider laws begins with 
an examination of the algebraic power-balance expressions 
in matrix-vector form. From these, we extract the exact 
nonlinear mapping between line complex-power flows and 
nodal complex-power injections in analytical closed form. The 
first step proceeds in the same vein as a, and it involves 
uncovering the contributions of nodal current injections to 
line-current flows. It emerges that this is a linear function of 
entries of the network admittance matrix and independent of 
the network’s voltage profile, i.e., bus-voltage magnitudes and 
phase angle differences. We refer to the sensitivities of line 
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current flows to nodal current injections as current injection 
sensitivity factors. Leveraging this first step, we then derive 
the impact of bus complex-power injections on line complex- 
power flows through algebraic manipulations of the power- 
flow expressions written in matrix-vector form. The resultant 
nonlinear expressions that constitute the power divider laws are 
functions of the voltage profile of the AC electrical network, 
as well as network topology (since they include the current 
injection sensitivity factors by construction). Additionally, we 
employ engineering insights on small angle differences, a 
uniform voltage magnitude profile, and the inductive nature 
of transmission networks to obtain simplified approximations 
of the power divider laws. For instance, one such simplification 
decouples the active- and reactive-power flows and injections. 
Under the most restrictive set of assumptions, power divider 
laws yield the ubiquitous DC power flow expressions jia. 

The remainder of this manuscript is organized as follows. 
Section im establishes mathematical notation and the power 
system model. In Section [Hll beginning with current injection 
sensitivity factors, we outline the derivation of the power 
divider laws. Approximations that emphasize decoupling of 
active- and reactive-power flows and injections as well as 
connections to the canonical DC power flow are provided in 
Section |IV] Section |V] highlights applications of the power 
divider laws to transmission-network allocation, loss alloca¬ 
tion, and feasible-injection identification. These applications 
are illustrated via numerical case studies in Section |Vll Fi¬ 
nally, concluding remarks and directions for future work are 
provided in Section IVIII 

II. Preliminaries 

In this section, we introduce relevant notation and describe 
the power system model used in the remainder of the paper. 

A. Notation 

The matrix transpose is denoted by (-j^, magnitude of a 
complex number by | • |; complex conjugate by (•)*, complex- 
conjugate transposition by (-j^, real and imaginary parts of a 
complex number or vector by Re{-} and Im{-}, respectively, 
and j := y/—!. A diagonal matrix formed with entries of 
the vector x is denoted by diag(a:); and diag(a:/y) forms a 
diagonal matrix with the m-th entry given by XmlVm, where 
Xm and yjn are the m-th entries of vectors x and y, respec¬ 
tively. The spaces of x 1 real- and complex-valued vectors 
are denoted by and C^, respectively, and denotes 
the A^-dimensional torus. The spaces of M x N real- and 
complex-valued matrices are denoted by and £^xn^ 

respectively. The entry in the m-th row and n-th column of 
the matrix X is denoted by [X]mn- The x 1 vectors with 
all ones and all zeros are denoted by and ©at; denotes 
a column vector of all zeros except with the m-th entry equal 
to 1; and Cmn ■= Sm — e„. For a vector x = [a;i,..., xat]"'", 
we denote cos(a:) := [cos{a;i),..., cos(xAr)]'^, sin(a:) := 
[sin(a;i),..., sin(a;Ar)]'^, and exp{a::) := • ■ •) Fi¬ 

nally, .— [ci, . . . , Cm —17 Cm-t-l 5 ■ ■ • 5 C^v] € R 
Corresponding to the vector x G , we denote a;"* := II^x; 
note that a;™ is recovered from the vector x by removing the 
m-th entry. 


B. Power System Network Model 

Consider a power system with N buses collected in the set 
J\f. Lines in the electrical network are represented by the set 
of edges E := {(m, n)} C Af x Af. Denote the bus admittance 
matrix hy Y = G -\-]B G , Transmission line (m, n) G 

E is modelled using the lumped-element Il-model with series 
admittance ymn = Vum = gmn+]bmn G C \ {0} and shunt 
admittance j&mn £ C \ {0} on both ends of the 

line. Entries of the admittance matrix are defined as 

f 2/m-l-E(m,fe)G£y”^fc> if m = n, 


[^]mn ■— \ ymn-, 

if (m, n) G E, 

(1) 

l 0: 

otherwise, 


where 



Vm — 9m 

■= 2/mm ^ ^ 2/mfcJ 

(2) 


feeA'm 

denotes the total shunt admittance connected to bus m with 
Afm C Af representing the set of neighbours of bus m and 
2 /mm G C any passive shunt elements connected to bus m. 
Notice that if the electrical network had no shunt elements, 
then y is a singular matrix with zero row and column 
sums. On the other hand, the inclusion of (even one) shunt 
admittance term in ([T]) intrinsically guarantees invertibility of 
Y by rendering it irreducibly diagonally dominant lITSll . 

Let V = [Vi,..., Vn]"^, where Vi = \Vi\Y9i G C represents 
the voltage phasor at bus i, and let / = [/i,..., where 

li G C denotes the phasor of the current injected into bus 
i. The bus-voltage magnitudes are collected in the vector 
\V\ = [|Vi|,..., |V)i|]'^ e R^q, and bus-voltage angles are 
collected in the vector 9 = [0i,... ,6*Ar]^ G T^. Kirchhoff’s 
current law for the buses in the power system can be compactly 
represented in matrix-vector form as 

I = YV. (3) 

Denote the vector of bus complex-power injections by 5" = 
[5'i,...,5'Ar]'^ = P + jQ, with P = [Pi,...,PAr]'^ and Q = 
[Qi, ■ ■ ■ ■ (By convention. Pi and Qi are positive for 

generation and negative for loads.) Then, bus complex-power 
injections can be compactly written as 

5 = diag(U)r. (4) 

The equation above is the complex-valued equivalent of the 
standard power flow equations expressed in a compact matrix- 
vector form, and generalized to include active- and reactive- 
power injections as well as voltage magnitudes and phase 
angles at all buses (including the artifactual slack bus). We 
will find the following auxiliary bus voltage phase angle- and 
active power-related variables useful; 

9^:=9mtN-9, 9"^:=Ul9^, P™ := P, (5) 

where 0 ™ is obtained by setting the system angle reference 

as 9m (the voltage phase angle at bus m). In 0 , and 

P™ denote the {N — 1)-dimensional vectors that result from 
removing the m-th entries from the W-dimensional vectors 0 ™ 
and P respectively. For example, with respect to the system 
illustrated in Fig. |2] where N = 3, relevant variables are 0 = 
[ 01 , 02 , 03 ]"^, P = [Pi, P 2 , Ps^^■ With the choice m = 2, it 
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follows that 0^ = {62 —91^^,92 — 07!^, = [6*2 —^i,^*2 —^3]"*", 

and P2 = [Pi,P3]T. 

III. Derivation of the Power Divider Laws 
In this section, we demonstrate how line current flows can 
be expressed as a linear function of bus curi'ent injection 
phasors based on entries of the network admittance matrix. 
Leveraging this relationship, we then derive expressions that 
uncover the contributions of bus active- and reactive-power 
injections to line active- and reactive-power flows. 

A. Line Current Flows 

We begin by expressing the curi'ent in line (m, n) S f as 

^{m,n) — ymn^Y^m ^n) IJm^Tn 

= {Umn^rnn ”1" ym^m) (6) 

Since the bus admittance matrix, Y, is invertible, from O, 
the bus voltages can be expressed as V = Y~^I. As a direct 
consequence, (|6ll can be written as 

I{m,n) ym^m) ^ ^ 

where K(^rn,n) G C^. The entries of referred to as 

the current injection sensitivity factors of line (m, n) with 
respect to the bus current injections. The current injection 
sensitivity factors in © delineate the exact effect of bus 
current injections on the curient in line (m, n). Moreover, they 
depend only on network parameters, and not the operating 
point, i.e., the voltage magnitudes and phase angles across the 
network do not influence ©. For subsequent developments, 
we will And it useful to decompose the current injection 
sensitivity factors into real and imaginary components as 

CI(m,n) “f j/5(m,n)- (^) 

Note that the expression in © is a generalization of the current 
divider law, which typically applies to the particular case of a 
single current source feeding a set of impedances connected 
in parallel ©. 

Remark 1 (Invertibility of Y). In ©, invertibility of Y is 
implicitly assumed. On the other hand, suppose Y is not 
invertible, i.e., there are no ground-connecting shunt elements. 
The current in line (to, n) is then given by 

^{m,n) — ymn{Ym ^n) — yrtin^rnnY- (^) 

In this case, the bus voltages satisfy 

V = Y^I+^lNtlV, (10) 

which follows from pre-multiplying both sides of © by the 
pseudoinverse of the admittance matrix, Y'^, and recognizing 
that the admittance matrix and its pseudoinverse are related 
by 113 

yyt = y^Y = diag(lAr) - 

Substituting for V from ( fTOl i in ©, we see that the current 
injection sensitivity factors in this case are given by 

Kjm,n) = (11) 

which follows from the fact that = ©^. ■ 


B. Line Active- and Reactive-power Flows 

Denote, by -f jQ(m,„), the complex power 

flowing across line {m,n). We can write 

^(™.n) = KnI(V)' (12) 

We substitute the current injection sensitivity factors from © 
(or (fTTl i for non-invertible admittance matrices) into (fT3 l. and 
obtain 

Sim,u) = (13) 

Eliminating I* from ([TtI i using ©, we get 

S(m,n) = (^iag S. (14) 

With the phasor form of the voltages, we can write GUI as 

= (iKnlei^™) «f^,„)diag (15) 

where we point out that with reference to the notation estab¬ 
lished in Section Hl-Al 

(r^)" ■ ■ ■ ’ ■ 

Since line power flows are often expressed as functions of 
angle differences, we And it useful to rewrite (ITsT i as 

S(m,n) = |Kn|Kf„,„)diag 

= | 14 ^|Kf^,„)diag 5', 

where 0™ = Om^N — 9 G T-^. To simplify the expression 
above, define: 

, /cos6l'"\ , /sin^^X 

Tjfl'm := diag I I +i diag I ) • 

Making use of the above and the decomposition of the current 
injection sensitivity factors in ®, we obtain 

S{m,n) = |14n| ^Ct(m,n) ~ + j4'm)(T’+ jQ) 

= |14m| + P(m,n)^ d" jQ) 

+ j|14n| ~ {P + jQ)- 

Then, we get the following power divider laws that indicate 
how bus active- and reactive-power injections map to the 
active- and reactive-power flows on line {m,n): 

P{m,n) lde{S'(m ji)} |T)ti| {f^{rn,n)P i (16) 

Q(m,n) — IlH{S'(m,n)} = \^m\ (^'^(m,n)Q d” ' (1^) 

In (GSIi-iGtI), U(^jn,n),V(m,n) G are given by 

tl(m,n) ^mO!(^rn,n) d“ ^mP(m,n) (1^) 

, /'cos6»'"\ , /sin6»™\ ^ 

= diag I^Tp J + diag j P(m.n), 





m,n) ■— 


= diag 


m,n) ^mP(‘ 

sin 9 


1^1 


m,n) 

a(m.n) - diag 


(19) 


cos 0" 

7^ 


/3(ri 
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We note that (fT6ll-(fT7li are the complex-power analogues of the 
current divider in O, namely they indicate how the complex- 
power injection at each bus i contributes to the complex-power 
flow on line (m, n). We refer to (fT6l l- (fT7l i as the power divider 
laws, since they specify how the active- and reactive-power 
flows on line (m, n) divide among active- and reactive-power 
injections at each bus i. 


The power divider laws in (l2Tt- (l24l i are then given by 

P{m,u) = (di&g p - diag , 

(27) 

Q(m,n) = (^diag Q + diag . 

(28) 


IV. Approximations, Decoupling, and Connections 
TO THE DC Power Flow 

In this section, we make use of a few practical observations 
and engineering insights regarding high-voltage transmission 
systems to present a suite of approximations to (HDl-linii 
that are conceivably applicable in different contexts. We also 
uncover the classical DC power-flow expressions under the 
most restrictive set of assumptions. 

A. Lossless Networks 


C. Unity Voltage Magnitude 

In the per-unit system, the numerical values of voltage 
magnitudes are typically such that \V\ ~ In- Thus, (l25]) - (l26l) 
simplify to 

~ 0^(m,n)5 (29) 

~ ding ) Cy{^m,n)- (^d) 

Accordingly, the approximate power divider expressions 
in are then given by 


The line resistance in transmission circuits is much lower 
than the line reactance lfT2ll . As a direct consequence, the con¬ 
ductance is much smaller than the susceptance. In other words, 
9mn << bmn in the series admittance of line (m, n) G £ and 
9m << bm in the shunt admittance at bus m G Af. Thus, we 
can approximate Ri j&„„, r; j6„, and Y r; jB. With 
these approximations in place, O simplifies to the following; 


^(m,n) ~ ^(m,n) ~ {bmnS-mn + bmCm) P ■ (20) 

Therefore, under the lossless assumption, (fT8l l-(fT9]l become 

" cos 0" 


't^(m,n) ~ ‘^mt^(m,n) diug 


1^1 

/ sin 0 ” 


t^(m,n)i (21) 


~ '^mt^(m,n) diug ^ j 0^(m,n)- 


( 22 ) 


The power divider expressions in (fT6l l- (fT7]) are then given by 

' cos 0 ” 


P(m,n) = (^diug ^ 


— diag 


Q{m,n) — |hm|ct(m,ra) ( diag 


diag 


1^1 

sin 

7^ 

cos 0 ^' 

77 

sin 

77 


p 

q) , (23) 


Q 

P] ■ (24) 


B. Small Angle Differences 

Under typical operating conditions, the phase angle differ¬ 
ences between any two transmission-system buses are small, 
i.e., 0i — 0k << 1, for all iff G Af. It follows that 
sm{0i — 0k) « 0i — 0k and cos{0i — 0k) ~ 1. Under these 
assumptions, EB-Elli further simplify to 






diag 

diag 







(25) 

(26) 


Pim,n) = ajm,n) (P “ ^iag (0"^) Q) , (31) 

Q(m.n) = ajm,n) (Q + diag (6''") P). (32) 


D. Decoupled Power Flow 

Since angle differences are expectedly small, it follows that 
eJ0'^Q^ « Pi and ej0^P^ « for all i G Af. With 
this assumption in place, the expressions in (l3Tl i- (l32l i further 
simplify to 

P(m,n) = ^{m,n)P^ (33) 

Q{m,n) = ^(m,n)Q- (34) 

The above expressions echo the widely accepted decoupling 
assumptions between active- and reactive-power in transmis¬ 
sion systems. Note that (l33]) - (l34l) are valid only when the 
entries of P and Q are of the same order of magnitude. For 
instance, if reactive-power injections are close to zero then the 
approximation in 0^ would no longer be valid. 

E. Connection to the DC Power Flow 

Here, we demonstrate that the DC power-flow equations can 
be derived from the active-power flow expression in (l3^ under 
the further approximation of neglecting shunt susceptance 
terms in a lossless transmission network. In order to uncover 
this relationship, a slack bus needs to be assigned (since 
the construction of the DC power flow equations requires a 
slack bus). Without loss of generality, set bus 1 as the slack 
bus. Accordingly, partition = [a(m,n).i, 

= eTo:(m,n) C R isolates the element of a(^m,n) that 
corresponds to the slack bus, and a(^rn,n) = 7ci(m,n) C 
|R^“t collects the remaining elements. Similarly, partition 
P = [Pi,P^], where P^ = IlJP. With these partitioned 
vectors, (l3^ becomes 

P(m^n) “f ^(^m,n)P 


(35) 
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where the second equality above results from substituting 
Pi = which is valid for lossless systems, since 

iJ^P = Pi+ lJf_iP'^ = 0. 

Next, we examine the structure of for a system 

in which losses and shunt elements are omitted. Neglecting 
shunt susceptance terms (i.e., setting bm = 0, Vm £ Af) 
in (Eil, we obtain ■ Note that the bus 

susceptance matrix B is not invertible since ground-connecting 
shunt elements are removed, so we resort to using the pseudo¬ 
inverse of B, namely B^ (see Remark[T]i. We can equivalently 
express n)^ ~ ^mnS^n- Partitioning out bus 1 (the slack 
bus) in a(m.n)> B, and emm we get that 


;rT 


biA bl 



1 


— ^mn ^m7 


(36) 

where B is constructed by removing the row and column 
corresponding to the slack bus from B. Focusing on the second 
element in (l36l l. we obtain 


^(m,n 




<^{m,n)B = b, 


mn^rni 


(37) 


Furthermore, since shunt elements are neglected, B satisfies 
BIm — ©jv, from which we get bi — —By 
substituting this into (IJTT i. it simplifies to 

'^{m,n) ~ Q!(m,rt),lljv-1 J ^ ~ bmn^mn 

^(m,n) ~ ~ bmnS-mnB ■ (38) 


where the second equality above results from recognizing 
that the reduced bus susceptance matrix B is invertible. 
Substituting (|3^ into (|35]l, we obtain 

P(7n,u) = bmn^nB-^P"- ( 39 ) 


The classical expression for the power flow on a line with a 
small-angle approximation is given by ED 

P{m,n) — ^n) — bmn^rnn^ ’ (40) 

where 0^ is as defined in (|5ll. Finally, by equating (l40l i 
with (|39]l, we get = —B~^P^, which are the DC power- 
flow expressions lEl. 


Example 1 (Approximations). Consider the 3-bus system with 
the one-line diagram shown in Fig. |2] Transmission lines are 
modelled with lumped parameters, where yi 2 = 1.3652 — 
jll.6041p.u., = jO.OSSp.u., 2/23 = 0.7598-j6.1168p.u., 

yil = j0.153p.u., yi 3 = 1.1677 - jl0.7426p.u., yg = 
j0.079p.u. In this system, generators are connected to buses 


TABLE I: Exact and approximate line active- and reactive- 
power flows for Example [T] All quantities are in p.u. 


Line 

Exact 

Lossless 

Small 

Unity 

DC 

Flow 



Angle 

Magnitude 

Power Flow 


iiD-ini 

(23)-(24) 

lUl-llD 

ISD-OU 

|40) 

Ul.2) 

0.0533 

0.0515 

0.0461 

0.0753 

0.0300 

P(2,3) 

0.844 

0.843 

0.843 

0.847 

0.800 

Pa.3i 

1.54 

1.55 

1.55 

1.52 

1.43 

Q(1,2) 

0.0821 

0.0894 

0.0880 

0.0965 

— 

Q(2,3) 

-0.0123 

-0.0061 

-0.0059 

-0.0051 

— 

<3(1,3) 

0.370 

0.363 

0.364 

0.356 

— 


1 and 2, injecting Pi = 1.5973 p.u. and P 2 = 0.7910 p.u., 
respectively. A PQ load is connected to bus 3, with P 3 = 
—2.35 p.u. and = —0.5 p.u. The voltage magnitudes at 
buses 1 and 2 are regulated at |Vi| = 1.04 p.u. and IV 2 I = 
1.025 p.u., respectively. Using the parameters listed above, we 
compute the power flow solution and then proceed to compute 
the line active- and reactive-power flows using the exact 
expressions in ([I 6 ]l and dnii, respectively. Subsequently, we 
approximate the line flows using the simplifying assumptions 
outlined in Section IIVI and results are collected in Table |I] ■ 


V. Applications 

In this section, we detail three applications for the power 
divider laws derived in Section |III] and their approximations 
described in Section HVl 


A. Transmission-network Allocation 


With the power divider laws, we can explicitly uncover the 
proportion of bus injections that contribute to line flows in ana¬ 
lytical closed form. In the literature, this is typically referred to 
as transmission-network allocation, and when electricity prices 
are factored in, as transmission-network cost allocation 13 
Erom (fTST l and (fTTI) . we can express the active- and reactive- 
power flows on line (m, n) as the following sum of 2N scalar 
terms; 


P( 


(m,n) 


Q (m,n) 


N 


N 


/ N N \ 

IK.I E t-i ^{7n,n)Qi T ^ ^ 6^ ^{m,ri)Pi 


\i=l 


i=l 


, (41) 

. (42) 


With (HTI) . for each bus i, we compute the contributions of its 
active- and reactive-power injection components to the active- 
power flow in line (m, n) respectively as 




Pi 


(m,n) 


(m,n) 


(43) 


Similarly, from (|42] |. we compute the contributions of the 
active- and reactive-power injection components at bus i to 
the reactive-power flow in line (m, n) respectively as 



Pi 


U(m,n)Qi 
Q {m,n) 


(44) 


Example 2 (Active-power Transmission Allocation). Consider 
the 3-bus system from Example [T] Without loss of generality, 
let us focus on the active-power flow on line (1,3), which 
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is P(i^ 3 ) = 1.5440 p.u. Computed using (l4?b . the bus active- 
power injection contributions to P(i^ 3 ) from buses 1 , 2 , and 
3 are 49.88%, 12.11%, and 39.19%, respectively. The bus 
reactive-power injection contributions are on the order of 
1%. These quantities can be directly used by power system 
operators to apportion the cost of transmitting active power 
on lines in the network to generators and loads. ■ 

B. Transmission-loss Allocation 

The expression in (fTbl) can be leveraged to uncover the 
contributions of bus active- and re active-power injections to 
active-power transmission losses on each line. In the literature, 
this general task is typically termed as transmission-loss 
allocation jSl. Denote the active-power loss on line {m,n) 
as T(m,ri)> and consider that T(m,n) can be expressed as 

\^{7n,n) ym^m\ ( 45 ) 

= Re {2/77171 (f(m,n) Unfi^va) 2/ml^)} ■ 

Recognizing from dUl that I(^rn,n) - VmVm = {Vm - Vn)ymn, 
and substituting this into ( l45l l. we obtain 

Lim,u) = Re {(^7 - Vr,)y*^„iV^ - V^)*} . ( 46 ) 

Further suppose the shunt elements are purely imaginary, i.e., 
ym = ibm, which is a common assumption in transmission 
line modelling, then (l46T l can be equivalently expressed as 


injections, specifically. Pi = 1.60 p.u., P 2 = 0.791 p.u., and 
P 3 = —2.35 p.u. ■ 

C. Line Active-power Flow-constrained Injections 

Here, we consider the problem of obtaining a set of bus 
active-power injections that best satisfy a set of desired 
transmission-line active-power flows. In order to achieve this 
goal, below, we formulate a convex optimization problem that 
utilizes the power divider approximation in Section IIV-DI 
Let V C £ denote the subset of D transmission lines for 
which active-power flows are specified and suppose D > N. 
Collect these D designated reference line flows, i.e., {PJ^ 

(m, n) G P into the column vector P-jj G and correspond¬ 
ing real components of current injection sensitivity factors, i.e., 
{a( 77 j 77 )}, (to, n) G V, into the matrix A G In order to 

obtain a set of bus active-power injections that best satisfies the 
desired line flows and simultaneously adheres to system power 
balance, we solve for P G R^ from the following linearly 
constrained least-squares optimization problem 01 : 

min \\AP-Pv\\^ 

pgr" 

s.t. llP = L, (49) 

where L denotes the total system loss. The unique closed-form 
solution to (|4^ is given by: 


L(m.u) = Re{(U7n - Vnjyl^n^Vm " K)* 
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'2Af^A In 

-1 

'2A^Pv 

+ y^yiy:a + ynyivvi , 

(47) 
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_ 0 


L 


(50) 


where the addition of Vmym^m + ^nyV^n^ affect 

the real part of the expression. Finally, with I(^rn,n) in ®, we 
simplify (l47l l as 

P{m,n) ^^rnl(^rn,n) “f P(m,n) P P{n,m)j 

where P(m,n) {P(n,m)) represents the active-power flow from 
bus TO (n) to n (to) along line (m,n). With the aid of (fThl l. 
we can write 



which indicates how bus active- and reactive-power injections 
contribute to the active-power line loss on the (to, n) line. 

Corresponding to each one of the cases in Section IIVI the 
expression for the transmission-network losses in (l48T l also 
simplifies correspondingly. In the particular case considered in 
Section HV-Dl the expression for the losses on the {m,n) line 
simplifies to = {c({m,n) + ct{n,m))P^ which resonates 

with prior art that has used current flows/injections as proxies 
for power flows/injections 0 . 

Example 3 (Transmission-line Losses). Consider, again, the 
3-bus system from Example [T] To demonstrate the concepts 
above, using (ESll, transmission-line losses are obtained for 
each line as T(i, 2 ) = 0.0003 p.u., T( 2 , 3 ) = 0.0140 p.u., 
and L(i 3 ) = 0.0240 p.u. The total system loss can be 
computed by summing up losses in all lines as 0.0383 p.u. 
We verify this quantity by summing up all bus active-power 


where A is the Lagrange multiplier associated with the equality 
constraint. Invertibility of the matrix in ( l50t can be ensured 
if and only if has independent columns ifTSl . Note 

that the total system loss L in (l49l l is not known a priori. One 
option is to set it to zero, as would approximately be the case 
in transmission networks with dominantly inductive lines. An 
alternative approximation is highlighted next. 

Neglecting shunt elements in (l45T l. we get that 
L{m,n) = l^(m.n)pR-e{ 2 /,j;i}. Recall, from Section IV-BI 

that current flows/injections can act as proxies for power 
flows/injections. Consequently, we can approximate 
L{m,n) ~ l*S'(m,n)pRe{y,j;( 7 }. Furthermore, suppose that 

the power factor on each line (to, n) is close to unity, 
then we can express the expected loss on line (to, n) as 
Hm,n) = l-P(m.n)PRe{y-j7}, where denotes the 

desired active-power flow on line {m,n). Similar to the 
construction of Px>, collect these predicted line losses, i.e., 
{'^(mn)J'’ G V, into the column vector G R^. 

With these variables in place, we reformulate (l49l l with 
L = l^Lu, resulting in the following unique closed-form 
solution: 

P 
A 

Example 4 (Constrained Injections). Consider the 3-bus sys¬ 
tem from Example [1] Using O and ([ 8 ]l, we get a(i 2 ) = 
[0.518,-0.233,0.249]^ a( 2 , 3 ) = [0.244,0.493,-0.0289]^ 
and a(i, 3 ) = [0.482,0.233, -0.249]"’’. Suppose the desired 
reference active-power flows on lines in the set V = 


2 A^A 

In 

— j. 

2 A^Pv 

. -“-jv 

0 






















7 



1.3 ^ 


0.8 

0.3 

- 0.2 


P 


o 

Oh 

-0.7 ^ 


-1.2 g 





10 


a 

o 


-0 S 


03 

CD 


■10 


(a) Active-power flow. (b) Reactive-power flow. 

Fig. 3: IEEE 14-bus system: Contribution of bus active- and reactive power (marked as circle and square, respectively) injections 
to active- and re active-power flow in line ( 6 ,12) (arrow indicates direction of flow). Each bus is associated with a coloured 
circle and a coloured square, with the colour representative of the bus’ contribution to the power flow in line ( 6 , 12 ). 


{(1,2), (2,3), (1,3)} are = 0.46 p.u., = 

0.67 p.u., and = 1.65 p.u.. We have thereby assembled 
all ingredients of the optimization problem in ( |49] |, where A = 
[a(i,2)ca(2,3):a(i,3)]'’". and Pp = [■P('i, 2 )c-^( 2 , 3 )’ 

Based on desired line active-power flows 
above, the expected line losses are collected in 
Lv = [-^( 1 . 2 )’-^( 2 , 3 )’-^( 1 , 3 )]'’"’ where = 0.0021 p.u., 

P }2 3 ) = 0.0090 p.u., and = 0.0272 p.u. Via dSTT l. we 

get that Pi =2.11 p.u., P 2 = 0.222 p.u., and P 3 = —2.29 p.u. 
With these bus active-power injections, the nonlinear power 
flow solution reveals that the actual line active-power 
flows are P(i. 2 ) = 0.468 p.u., P( 2 , 3 ) = 0.688 p.u., and 
^*( 1 . 3 ) = 1.64 p.u. The 2-norm of the deviation of the 
actual line active-power flows from the desired ones is 
0.0218 p.u. The actual system total loss is 0.0384 p.u., 
which is approximately equal to the expected value of 
= 0.0383 p.u. 

Next, neglecting losses in (ISOl l with L = 0, we get Pi = 
2.11 p.u., P2 = 0.208 p.u., and P3 = —2.32 p.u. With these 
updated bus active-power injections, we compute the nonlinear 
power flow solution and And that the actual line active-power 
flows are P(i.2) = 0.486 p.u., P(2,3) = 0.692 p.u., and 
P(i, 3 ) — 1-66 p.u. The 2-norm of the deviation of the actual 
line active-power flows from the desired ones is 0.0360 p.u. 
Higher errors are attributable to neglecting losses. ■ 

VI. Case Studies 

In this section, we illustrate concepts presented in Sec¬ 
tions HIIHV] using the IEEE 14-bus system. The simulation tool 
MATPOWER iflbl is used throughout to compute all power- 
flow solutions. 

A. Transmission-network Allocation 

In this case study, as depicted in Eig. [3 we visualize the 
proportion of bus active- and reactive-power injections that 
contribute to active- and reactive-power flows in line ( 6 , 12 ). 
In Eig. 1^ the circle at each bus i represents the percent 
contribution of its active-power component, i.e., the fth term 
in IVeju^ 12 )-^’ active-power flow in line ( 6 , 12 ), and 



Fig. 4: IEEE 14-bus system: Contribution of bus active and 
reactive power (marked as circle and square, respectively) 
injections to active-power loss on line ( 6 , 12 ). 

the square represents the percent contribution of its reactive- 
power component, i.e., the i-th term in —\Vq\vJq-^ 2 )Q^ both 
computed via (l43l l. Similarly, in Eig. the circles and 
squares indicate the bus active- and reactive-power injection 
contributions, respectively, to the reactive-power flow on line 
( 6 ,12), computed via (l44ll . As described in Section IIV-DI 
we note that, indeed, there exists strong coupling between 
the line active-power flow and bus active-power injections. 
On the other hand, the coupling is not evident between the 
line reactive-power flow and bus reactive-power injections. As 
mentioned in Section IIV-DI this is because the power factors 
of complex-power injections at most buses are close to unity. 

B. Transmission-loss Allocation 

In this case study, we uncover the proportion of bus active- 
and reactive-power injections that contribute to active-power 
loss in line (m, n) G E. We compute transmission losses for 
each line (m, n) in the IEEE 14-bus system via the expression 
in dull. Similar to (HTt . we can rewrite (l48T l as the sum of 
2N terms and express, for each bus i, the contribution of its 
active-power injection to the loss on line (m, n) as 

ef {\Vm\u (m,n) + \Vn\U( n^m ))J\ 


(52) 

































Fig. 5: IEEE 14-bus system; Histogram of 5000 line active- 
power flow-constrained injection simulations. The horizontal 
axis represents the 2-norm of the error between the actual line 
active-power flows and the desired ones, while the vertical 
axis represents the frequency with which a particular deviation 
range occurs. 

Analogously, for each bus i, the contribution of its reactive- 
power injection to loss on line (to, n) can be expressed as 

ef {\Vm\V( m,n) + Qi 

j • 

In Eig. |4] we plot the bus active- and reactive-power injection 
contributions to the active-power loss on line (6,12). Note 
that, unlike for the line active-power flow, there is no clear 
decoupling between active-power loss in line (to, n) and bus 
reactive-power injections. Eor example, 27.4% of the loss in 
line (6,12) is contributed from the active-power injection at 
bus 14 and —16.8% is from the reactive-power injection at bus 
13. The negative sign is indicative of injection directionality 
(i.e., generation versus load). 

C. Line Active-power Flow-constrained Injections 

In this case study, using the IEEE 14-bus test system, we 
illustrate the application of the power divider laws to obtain a 
set of active-power injections that best satisfy a set of desired 
transmission-line active-power flows. To simulate these desired 
reference line flow values, for each line (to, n) S E, we apply 
random perturbations to the base-case transmission-line flows: 

P(m,n) = P{m,n) (1 + O') i (54) 

where cr is a random variable uniformly distributed between 
— 1 and -t-1. With these randomly generated desired flow 
values, both lossy and lossless assumptions are applied and 
corresponding solutions are obtained via (EB and (fSOl l with 
L — 0, respectively. Eor each set of resulting bus active-power 
injections, the nonlinear power flow solution is computed to 
obtain the actual line active-power flows. These are then com¬ 
pared to the desired reference values, and the 2-norms of errors 
between the actual flows and the desired ones are recorded. In 
total, we run 5000 simulation cases and plot the histogram of 
errors in Eig. |5] We note that the lossless approximation leads 
to line flows that are closer to the desired reference values. 


This is likely due to the fact that the magnitude of the complex- 
power flows cannot be approximated by their real components 
only, i.e., the power factors of actual flows are not sufficiently 
close to unity. Einally, we mention that it is quite likely that 
the randomly generated desired set of line active-power flows 
is, in fact, infeasible, and hence may lead to actual line flows 
with large deviations from the desired ones. 

VII. Concluding Remarks and Directions eor 
Euture Work 

In this paper, we derive the so-called power divider laws, 
which indicate how active- and reactive-power flows on a 
transmission line are divided among bus active- and reactive- 
power injections. A suite of approximations are then derived 
based on simplifying assumptions commonly applied to the 
analysis of transmission networks. We present three applica¬ 
tions of the power divider laws: transmission-network allo¬ 
cation, transmission-loss allocation, and line flow-constrained 
injection estimation. Compelling avenues for future work are 
to leverage power divider laws for applications such as spot 
pricing, transmission-services pricing, and visualization. 
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